NON-PERTURBATIVE INVESTIGATION OF THE FERMION-HIGGS SECTOR 

OF THE STANDARD MODEL 
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ABSTRACT 

We present results for the renormalized quartic self-coupling Xr and the renormalized Yukawa coupling ya in a fermion- 
Higgs model with two SU(2) doublets, indicating that with the standard lattice regulator these couplings cannot become 
very strong. 



The Lattice Model 

It is an important issue to investigate within a non- 
perturbative regularization scheme whether the quar- 
tic self-coupling and Yukawa coupling of the fermion- 
Higgs sector of the Standard model (StM) may become 
strong when increasing the bare couplings to very large 
values while keeping the cut-off roughly constant. For 
this it is desirable to construct a lattice fermion-Higgs 
model with a fermion content which is as close as possi- 
ble to that of the StM. This is a non-trivial issue since 
a naive transcription of the continuum lagrangian with 
one SU(2) doublet (with gauge interactions switched 
off for simplicity) leads to the large number of 16 dou- 
blets on the lattice because of the species doubling phe- 
nomenon. Two strategies have been worked out so far 
to reduce this large number of mass-degenerate SU(2) 
doublets to one: The mirror fermion model [1] and the 
reduced staggered fermion model [2]. In this contribu- 
tion we shall only focus on the second method. 
Reduced staggered fermions are represented on the lat- 
tice by the one-component real Grassmann field \ which 
describes in the continuum limit two Dirac flavors. It is 
not immediately clear how to access these flavor degrees 
of freedom since the staggered spin and flavor degrees 
of freedom are spread out over the lattice. By using 
however the method outlined in ref. [2] it is possible to 
couple these two staggered flavors as an isospin-doublct 
to the Higgs field. For more details about the stag- 
gered method we refer the reader to refs. [2,3,4]. The 
final form of the one-doublet action in terms of the stag- 
gered field x reads 
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where Tp^ — y-r £ b ipn, x -b is the average of the scalar 
field over a lattice hypercube and n^ x = (— l) X1 ^ ^f- 1 , 
(pre = (— l) x f*+iH Vxi are the usual staggered sign fac- 
tors. The complete action is given by S = Sf + Sh, 



where S H = £ x [2«;£ 

^Pa X ^Pa, X -\-fi ^Pol X ^Pol X 

— ^{fa X fa X — l) 2 ] is the pure scalar field action. The ac- 
tion S is invariant under the so-called staggered fermion 
(SF) symmetry group, which contains a discrete set of 
0(4) flavor transformations. 

S is however not invariant under the full 0(4) flavor 
group: There are two operators with dimension four 
which are generated by the quantum fluctuations and 
which are invariant under the SF symmetry group, but 
break 0(4), 0« = £ « (2) = 



j Xl iVr^,x+p. 



Pij.cc) ■ I n order to recover the full 0(4) symmetry one 
has in principle to add these operators as counterterms 
to the action S — > S + eqO^ + SoO^ 2 ' and tune their 
coefficients £o and Sq as a function of the bare parame- 
ters such that the 0(4) invariance gets restored in the 
scaling region. We showed however in ref. [4] by using 
renormalized perturbation theory that the deviations in 
the renormalized couplings due to the 0(4) symmetry 
breaking (SB) are very small. The results presented in 
the following section have been obtained without adding 
these counter terms to the action. 
Since we are interested in the largest possible renor- 
malized couplings we have fixed in the numerical simu- 
lation X — oo. For the use of the Hybrid Monte Carlo 
algorithm it is furthermore necessary to use two mass- 
degenerate doublets. The n-y phase diagram is shown 
in fig. 1. There are four different phases, a paramag- 
netic (PM), a broken or ferromagnetic (FM), an an- 
tiferromagnetic (AM) and a ferrimagnetic (FI) phase. 
The various symbols mark the points in the FM phase 
where we carried out numerical simulations on lattices 
ranging in size from 6 3 24 to 16 3 24. 



Results of the Simulation 

The renormalized couplings have been computed from 
the usual tree level relations yn = mp/vn and = 
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bound 



Figure 1: Upper figure: Phase diagram at A = oo. 
Lower figure: m a /vn as a function of mp/oji. 

m 2 /2vft, where the renormalized field expectation value 
is defined as vr — vj\fZ^. Here v is the unrenormal- 
ized scalar field expectation value and Z v is the wave- 
function rcnormalization constant of the Goldstone prop- 
agator. For the determination of the quantities mp, 
ra a and Z^ we have measured the fermion-, er-particle 
and Goldstone-propagators in momentum space. The 
fermion propagator could be well described for all mo- 
menta by a one pole ansatz, which is characteristic for 
weakly interacting fcrmions. The numerical results for 
the cr-particle and Goldstone-propagators have been fit- 
ted to the ansatz G~^(p) = (p 2 + ml 7r + ^ sub (p))/Z a ^, 
where T, su b is the subtracted one fermion loop self- 
energy and p 2 is the lattice momentum squared. We 
kept m T as a fit-parameter, since the 0(4) SB effects 
give rise to a small non-zero value for m^, which how- 
ever is much below the other masses. The inclusion of 
the self-energy turns out to be essential if one wants to 
obtain reliable results for m^, and Z„^ also on smaller 
volumes. The fact that the one-loop ansatz is sufficient 
to describe the numerical results perfectly over a large 
momentum interval indicates already that the renor- 
malized couplings arc small. 

As a next step we have to extrapolate the finite volume 
results for yn and Xr to the infinite volume. We carried 
out the simulations on lattices of size L 3 24 with L rang- 



ing from 6 to 16. For the extrapolation to the infinite 
volume we use a 1/L 2 ansatz which is expected to hold, 
if the spectrum contains massless Goldstone bosons. In 
our case the Goldstone bosons in the infinite volume 
limit acquire a small mass due to the 0(4) SB. The fact 
that we have not observed significant deviations from 
the 1/L 2 behavior gives further evidence that the SB 
effects are small. 

In the lower graph of fig. 1 we display the infinite volume 
results for the ratios m a /vR = \ / 2Xr and mp/vR = yR. 
The symbols in the upper and lower diagrams of fig. 1 
match, so that one can see where in the phase diagram 
the results for the ratios have been obtained. It can be 
seen that the numerical values for neither ratio change 
when lowering n beyond k = 0. All points have roughly 
the same value of the cut-off in units of the scalar field 
expectation value vr rs 0.15 — 0.25. The arrows in fig. 1 
mark the tree level unitarity bounds for Xr and yR. The 
graph shows that the points obtained in the regions (II) 
and (III) of the phase diagram (see fig. 1) are still very 
close to these values, which indicates that the renor- 
malized couplings are not very strong. The solid line 
encloses the allowed regions obtained by integrating nu- 
merically the one-loop /3-functions from the cut-off scale 
down to the renormalization scale and identifying the 
couplings at these scales with the bare and the renor- 
malized couplings. The cut-off was adjusted such that 
the agreement with the numerical data is best. It is 
remarkable that the shape is in reasonable agreement 
with our data. Fig. 1 shows that the Yukawa inter- 
action gives only a modest increase in Xr. All in all 
we conclude that the renormalized quartic and Yukawa 
couplings are in accordance with triviality and that they 
cannot be strong, unless the cut-off is unacceptably low. 
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